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In previous papers I have argued that the fusion rules hypothesis, which was originally introduced by L'vov 
and Procaccia in the context of the problem of three-dimensional turbulence, can be used to gain a deeper 
insight in understanding the enstrophy cascade and inverse energy cascade of two-dimensional turbulence. In 
the present paper we show that the fusion rules hypothesis, combined with non-perturbative locality, itself a 
consequence of the fusion rules hypothesis, dictates the location of the boundary separating the inertial range 
from the dissipation range. In so doing, the hypothesis that there may be an anomalous enstrophy sink at small 
scales and an anomalous energy sink at large scales emerges as a consequence of the fusion rules hypothesis. 
More broadly, we illustrate the significance of viewing inertial ranges as multi-dimensional regions where the 
fully unfused generalized structure functions of the velocity field are self-similar, by considering, in this paper, 
the simplified projection of such regions in a two-dimensional space, involving a small scale r and a large scale 
R, which we call, in this paper, the (r,^)-plane. We see, for example, that the logarithmic correction in the 
enstrophy cascade, under standard molecular dissipation, plays an essential role in inflating the inertial range 
in the {r,R) plane to ensure the possibility of local interactions. We have also seen that increasingly higher 
orders of hyperdiffusion at large scales or hypodiffusion at small scales make the predicted sink anomalies more 
resilient to possible violations of the fusion rules hypothesis. 



I. INTRODUCTION 

A very surprising property of Navier-Stokes turbulence in 
three-dimensions is that it has an anomalous energy sink. This 
means that in the forced-dissipative case, with energy injected 
by random forcing at a constant rate the rate of energy dis- 
sipation e caused by the viscous term vV^Ua of the Navier- 
Stokes equation will eventually be equalized with even 
under the limit v — > 0+ where the viscosity v vanishes. This is 
surprising because one might expect that in the limit v 0+, 
the Navier-Stokes equation should lose its capability to dissi- 
pate energy. This does not occur because the vanishing vis- 
cosity is being compensated for by steeper velocity gradients. 

The existence of an anomalous energy sink in three- 
dimensional turbulence has been well-established by experi- 
ments and numerical simulations fl-Ht]. However, there is no 
proof, directly from first principles. An even more interest- 
ing situation exists for the case of two-dimensional turbulence 
(|4j-l6|]. In two dimensions, we have two conserved quadratic 
quantities, energy and enstrophy. Under the forced-dissipative 
scenario, again, with two dissipation terms, (a diffusion or 
hyperdiffusion term acting at small scales and a hypodiffu- 
sion term acting at large scales) we expect that most of the 
injected enstrophy will form an enstrophy cascade and will be 
dissipated at small scales and most of the injected energy will 
form an inverse energy cascade and will be dissipated at large 
scales. An open question is whether or not two-dimensional 
turbulence has an anomalous enstrophy sink at small scales 
and an anomalous energy sink at large scales, such that each 
sink does not contaminate its coiTesponding cascade. 

CuiTently, there is a lot of interest in understanding this 
question of anomalous sinks in terms of the singular solutions 
of the Euler equations |7]. However, it is also possible to ex- 
amine the same problem from the standpoint of the Navier- 
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Stokes equations. Consider, as an example, the case of the 
energy cascade of three-dimensional turbulence. Keeping the 
rate of energy injection Em constant, as we decrease the vis- 
cosity V, the energy cascade adjusts by moving the dissipation 
length scale 7], also known as the Kolmogorov microscale, fur- 
ther out into larger wavenumbers (i.e. smaller length scales). 
This intensifies the Laplacian V^m^ thus compensating for the 
decrease in the viscosity v at the dissipation term vV^mq.. In 
textbooks (e.g. in Ref. |j^), the location of the dissipation 
scale 77 is consequently derived by first assuming the existence 
of an anomalous energy sink and the existence of an energy 
cascade, and then deducing where the dissipation scale must 
be placed to dissipate the injected amount of energy. Thus, we 
obtain the well-known estimate r\ ^ (v^/e)^''^. 

What if we could turn this argument on its head? What if 
we could derive from first principles the location of the dis- 
sipation scale 7] without assuming a priori the existence of 
an anomalous energy sink, and then use that to prove that the 
energy sink is indeed anomalous? A rigorous form of this 
argument would have to establish from first principles: (a) 
that an energy cascade exists with energy spectrum scaling 
between k^^ and k'^ ' ; (b) the location of the dissipation scale 
itself. From (a) and (b) we may then conclude that there is an 
anomalous energy sink. This is obviously not an easy task; 
showing (a) is essentially almost the same thing as solving the 
problem of turbulence itself from first principles! Neverthe- 
less, we would still be interested in a mathematical argument 
that tolerates the making of some hypotheses, but then estab- 
lishes successfully both parts (a) and (b) without assuming, as 
Frisch |8]) did, a priori, the existence of an anomalous energy 
sink. 

Obviously, the problem of investigating the existence of an 
anomalous enstrophy and energy sink for the case of two- 
dimensional turbulence can be approached in a similar way. 
However, this task is complicated by the apparent lack of 
robustness of the cascades of two-dimensional turbulence, 
briefly discussed below, for which there is not currently a 
widely-accepted explanation @, I9l- [l5ll . For example, it is 
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well-known, from the work of Tran and Bowman ifTsUlsll . 
that the enstrophy cascade can fail to develop in the absence 
of a sufficiently strong large-scale sink. The large-scale sink 
is needed to dissipate the excess injected energy that cannot 
be disposed of by the small-scale sink. In the absence of a 
sufficiently strong large-scale sink, both energy and enstro- 
phy will pile up in the vicinity of the forcing wavenumber 
where they will eventually be dissipated when the pile-up be- 
comes sufficiently high. Likewise, for the inverse energy cas- 
cade, as Danilov and Gurarie il0l - [l2ll have shown, the inverse 
energy cascade can be "buried" under coherent vortices that 
hoard most of the energy of the system, resulting in a k^^ 
contribution to the energy spectrum that dominates the k^^^^ 
contribution of the still-existing inverse energy cascade. In 
a recent paper 1 16], 1 proposed that these coherent structures 
originate from a self-amplification of the sweeping interac- 
tions part of the Navier-Stokes non-linearity. So far as we 
know, the inverse energy cascade coexists with these coher- 
ent structures, and can be recovered by artificially eliminating 
them |9j A theoretical explanation for this recover- 

ability of the inverse energy cascade was given in previous 
papers SIHIll. 

Contrast this behavior with the case of the energy cascade 
of three-dimensional turbulence. While it is also possible to 
have a failure to develop an energy cascade, in a forcing- 
dissipation configuration whereby there is insufficient sepa- 
ration between the forcing scale and the dissipation scale, this 
can always be rectified by either increasing the rate of energy 
injection, or by decreasing the viscosity. In two-dimensional 
turbulence, a far more delicate tuning between forcing and the 
two sinks is required to recover a steady-state inverse energy 
cascade or a direct enstrophy cascade. We are thus faced with 
the additional challenge of accounting for the conditions that 
are required for the existence of the cascades, in addition to 
addressing the existence of anomalous sinks, when these con- 
ditions are satisfied. 

Further confusion arises from the fact that under the Tran- 
Bowman scenario ifTsUTsll . we can have a trivial type of an 
anomalous energy sink simply by having energy pile up at 
the forcing range until the dissipation terms become strong 
enough to balance the rate of energy injection. In fact, it can 
be shown that, at steady state, for finite viscosity, the total 
energy always has an upper bound |20]. It follows that two- 
dimensional turbulence always has an anomalous energy sink, 
even if it is merely of this trivial type. There is no known proof 
that a similar upper bound exists for the total enstrophy. As 
Eyink |20] noted, if one hypothesizes the existence of an up- 
per bound for the total enstrophy, one may then predict that the 
placement of the dissipation scales gives non-trivial anoma- 
lous sinks. Thus, when we discuss the existence of anomalous 
sinks in the context of two-dimensional turbulence, we should 
qualify that our main interest is in the existence of a non-trivial 
anomalous energy sink that dissipates the energy in a dissipa- 
tion range, located far away from the forcing range, so that 
the formation of an inverse energy cascade can be facilitated. 
We are, likewise, interested in the existence of a non-trivial 
anomalous enstrophy sink, again located far away from the 
forcing range, that can facilitate an enstrophy cascade. 



In earlier papers ifTsi [T9i I21I1 . I proposed that the local- 
ity and stability of the cascades of two-dimensional turbu- 
lence can be understood by employing a mathematical frame- 
work that consists of the generalized balance equations gov- 
erning the generalized structure functions F„ and the fusion 
rules hypothesis. This framework was originally introduced 
by L'vov and Procaccia in the context of studying the energy 
cascade of three-dimensional turbulence ||22| - [33|] . I proposed 
that the method underlying this approach could be general- 
ized to tackle the open questions that plague two-dimensional 
turbulence. In this framework, the generalized structure func- 
tions F„ are defined as products of velocity differences where 
each velocity difference is evaluated at two points in space dis- 
tinct from any other point-pair associated with the rest of the 
velocity differences. The usual standard structure functions, 
on the other hand, use the same point-pair for every veloc- 
ity difference. The fusion rules govern the scaling of these 
generalized structure functions when some but not all of the 
velocity point separations approach each other, while still re- 
maining within the inertial range. 

The structure of the argument within this mathematical 
framework, in broad strokes, runs as follows: We begin with 
the hypothesis that there is a region of scales in which the gen- 
eralized structure functions satisfy incremental homogeneity, 
incremental isotropy, and a reasonably weak hypothesis of 
self-similarity, and the requirement that these symmetries be 
universal. The fusion rules are derived from the universality 
hypothesis. Then, we proceed, with input from the govern- 
ing Navier-Stokes equations, to show that: (1) within this re- 
gion of scales the nonlinear interactions are local; (2) the self- 
similar scaling is indeed stable with respect to perturbations 
on the stochastic forcing term for the inverse energy cascade, 
and marginally stable for the downscale enstrophy cascade; 
(3) the dissipation scales are positioned as is necessary to pro- 
vide for anomalous sinks that can dissipate the injected en- 
ergy and enstrophy at wavenumbers far away from the forcing 
range. Parts (1) and (2) were shown in a previous paper ll2Tll . 
Part (3) is the subject of the present paper The entire argu- 
ment has been carefully summarized in the conclusion of the 
present paper 

The main thrust of the argument is to derive the location 
of the dissipation scales both at the downscale range and the 
upscale range. We then use that to conclude that the scaling 
dependence of the upscale energy dissipation e,> on the large- 
scale dissipation viscosity coefficient j3 is given by 

£,-,- ^ j3i-(&+2'n)/(f2-fe,l+2m)^ (j) 

Likewise the dependence of the downscale enstrophy dissipa- 
tion rate 77,,,, on the viscosity coefficient v is given by 

r?„v - vi-(f2-2(^+i))/(&,i-2(^+i))[in(£g/^^^,)]«3-i. (2) 

Here io is the forcing scale and A,,,, is the enstrophy dissipa- 
tion scale. Furthermore, 03 is a scaling exponent associated 
with the logarithmic scaling of the third-order vorticity struc- 
ture function, ^2 is the scaling exponent associated with the 
second-order velocity generalized structure function, £,2,1 is 
the fusion scaling exponent of that same structure function, K 
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is the order of the dissipation operator at small scales, and m 
the order of the dissipation operator at large scales. Under the 
fusion rules hypothesis we have ^2,1 = C2 for the downscale 
range, and ^2,1 = for the upscale range. For these values, 
the dissipation rates Tjj,,, and e,> become independent of the 
viscosities v and j3 . Thus, under the fusion rules hypothesis 
we have non- trivial anomalous sinks. We also see that increas- 
ing K and m leads to asymptotic independence of rjuv and e„ 
from V and /3 even when the fusion scaling exponents deviate 
from the prescribed values. 

Combined with the results of my previous paper |2ll], we 
are beginning to see a big picture in which the fusion rules 
hypothesis operates as a unifying nexus that can subsume 3 
distinct assumptions that everyone makes about turbulence 
cascades: locality, stability, and the existence of anomalous 
sinks. This unification provides us with two opportunities. 
First, it is possible to investigate the validity of the fusion 
rules hypothesis with numerical simulations, which are eas- 
ier to do in two dimensions. Second, and more important, 
the validity of the fusion rules can be investigated further on 
theoretical grounds, by generalizing similar arguments prof- 
fered for the problem of three-dimensional turbulence ||22| - 
I25I1 and the passive scalar problem [34]. It is worth not- 
ing that experiments have corroborated the fusion rules both 
for three-dimensional turbulence i35l - l40ll and for the passive 
scalar problem ll4lll42ll. 

This paper is organized as follows. Section 2 reviews the 
framework of the generalized balance equations and the fu- 
sion rules hypothesis. In Section 3 we show how the law gov- 
erning the location of dissipation scales follows as a conse- 
quence of the fusion rules hypothesis. Underlying this argu- 
ment is a geometrical conception of either inertial range as a 
two-dimensional region separated by the corresponding dis- 
sipative region, by a curve whose shape is deduced from the 
fusion rules hypothesis and locality (itself a consequence of 
the fusion rules hypothesis). The curvature of that line pro- 
vides us with an admissibility condition and also leads to the 
observable dissipation scale which is relevant to the anoma- 
lous sink question. In Section 4 we apply the method of Sec- 
tion 3 for the case of the enstrophy cascade, and in Section 5 
we consider the case of the inverse energy cascade. Having 
derived the laws governing the observable dissipation scales 
from the fusion rules for both cascades, in Section 6 we turn 
to the question of anomalous sinks. Conclusions and discus- 
sion, summarizing the logical structure of the argument as a 
whole, are given in Section 7. Technical matters are taken up 
in the appendices. 



II. SUMMARY OF PRIOR RESULTS 

We begin with a brief overview of the generalized bal- 
ance equations and the fusion rules hypothesis. The gener- 
alized balance equations were originally derived by L'vov and 
Procaccia ll28ll and, combined with the fusion rules hypoth- 
esis, they are the foundation of the non-perturbative L'vov- 
Procaccia theory of three-dimensional turbulence |26- 28, 30- 
[32I1 . and also the foundation for a corresponding investigation 



of the cascades of two-dimensional turbulence HIM in. 



A. The balance equations 

We begin with the Navier-Stokes equation for the velocity 
field in two dimensions: 



du 



(3) 



Here, Ua is the Eulerian velocity field, T^p is projection 
operator = d^p — dadj^V^^, fa is the forcing term, and 
D is the dissipation operator given by 



Dee (-l)''+'vV2'^ + (-l)'"+')3V 



m+ 1 o \j—2m 



(4) 



Here the integers K and m describe the order of the dissipation 
mechanisms, and the numerical coefficients v^- and j3 are the 
corresponding viscosities. D is the overall dissipation opera- 
tor The first term is the small-scale sink, and the second term 
is the large-scale sink. The case K" = 1 corresponds to stan- 
dard molecular viscosity, and the case m = corresponds to 
Ekman damping. 

From the Navier-Stokes equations we derive the exact sta- 
tistical theory of velocity differences. Let Wa(x,x',f) be the 
Eulerian velocity differences, defined as: 



Wa(x,x',f) = Ma(x,f) -Ma(x',f) 



(5) 



To write equations concisely, we introduce the following no- 
tation to represent aggregates of position vectors 



X-(x,x'), 
{X}„ = {Xi,X2, . . . ,x„}, 

{X}* = {Xi,...,X|(._i,X 



.,x„}. 



(6) 
(7) 
(8) 



Below, the notation ||{X}„|| ^ R means that all point to 
point distances in the geometry of velocity differences {X}„ 
have the same order of magnitude R. Similarly, the notation 
||{X}„|| <C ||{Y}„|| means that all the point to point distances 
in {Y}„ are much larger than all the point to point distances 
in {X}„. 

The eulerian one-time fully unfused correlation tensors are 
formed by multiplying n velocity differences Wa(x,x',f) eval- 
uated at 2n distinct points 



F„({X},„f) = 



k=l 



(9) 



Differentiating F„ with respect to t and applying the Navier- 
Stokes equations yields equations of the form 



dF 

-^+D„ = vJ„+liH„ + Q,„ 



(10) 



where D„ is the combined contribution of the pressure and 
the nonlinear term, Q„ is the contribution of the forcing term, 
J„ accounts for diffusion or hyperdiffusion, and H„ accounts 
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for large-scale dissipation. These are the generalized balance 
equations. The dissipation terms are given by 



A:=l 



(11) 



k=l 



where V^^ differentiates with respect to x^, and similarly with 

-2, 



V^,'^, V„,^'" and The forcing contribution is given by 



e„({x}„,f) = £e,„({x}^,Xi,f), 



k=l 



ete({x}„-i,Y,f) = 



A:=l 



(12) 



<Pj3(Y,f) 



where 

(Pa(X,f)=/a(x,f)-/a(x',f)- 

The non-linear term D„ can be rewritten as 

D„({X}„,f) = f D,„({X}„,f)+/„({X}„,f), 

<r=l 

with /„ representing the sweeping interactions given by 

/«({x}„,0 



(13) 



(14) 



= E('5i3,x, + ^i3,xv)(li/3({X}„,f) 



i:=l 



k=l 



where IXftdXln,?) is defined as 



1 # 



lt/j({X}„,f) = ^ L (Ma(Xi,0 + «a(x'i,0) , 



(15) 



(16) 



and Ditn representing the local non-linear interactions given, 
in general form, via a linear integrodifferential operator 0, as 
follows: 

£>^«({x},„0 

= //flfYifl'Y2 0(Xt,Yi,Y2)F„+i({X}f„,Yi,Y2,f). (17) 



The detailed form of the local nonlinear term D^^,, was given 
in our previous paper ll2lll . 



B. Self -similar scaling and the fusion rules 

For a stationary problem, the generaUzed balance equations 
can be rewritten as 



As we have explained in our previous paper 12111 , to con- 
siderable detail, the downscale or upscale inertial range of 
two-dimensional turbulence can be identified with a region 
3n C M?" in which for {X}„ e 3„, the contributions of the non- 
linear term 0„F„+i dominate all the other terms of the balance 
equations. The region 3n is circumscribed on one side by the 
forcing term Q„ and on the other side by the corresponding 
dissipation terms vJ„ and j3//„. Meanwhile, we assume, for 
the time being without further proof, that the sweeping term 
/„ also remains negligible with respect to 0„f„+i in the region 
3n- A preliminary study of the sweeping term was given in a 
previous paper lUql . 

Within the region 3n we expect that F„ will be self-similar 
according to the following scaling law: 



F„{X{X},„t)=X^"F„{{X}„,t 



(19) 



This scaling law is expected to hold when {X}„ G 3n, and the 
scaling exponents ^„ characterize the corresponding inertial 
range, and will obviously differ between the upscale range and 
the downscale range. Within the same region 3n it is expected 
that F„ will be incrementally homogeneous and incrementally 
isotropic. 

We also introduced [21] the hypothesis that both cascades 
of two-dimensional turbulence, when they exist, will satisfy 
another group of self-similarity laws called the fusion rules. 
We would now like to briefly summarize the content of these 
rules. Consider a geometry of velocity differences {x}„ such 
that all point to point distances have order of magnitude 1 , and 
define 

F,l''\r,R) = F„{r{xk}l„R{^,}Up+i)- (20) 

The function fi'''(r,/?) reflects the case where p velocity dif- 
ferences have separations with order of magnitude r, and 
n — p velocity differences have separations with order of 
magnitude R. The case of interest is when the evaluation 
{r{xi(}1^i,R{xk}'^^^j_^) is within the inertial range 3„ and 

r -^R. The fusion rules give the scaling properties of Fn''\ 

We distinguish between two cases. For a direct cascade, 
such as the energy cascade of three-dimensional turbulence 
and the enstrophy cascade of two-dimensional turbulence, the 
fusion rules are given by 



F^^ (Air, Aj/?) = Af '-"aI" ^'•■"i^^''^ 



(21) 



0„Fn+ i+I„ = vJ„ + ^H„ + Q„ . 



(18) 



for 2 < p < n — 2 with ^„ p = ^p. In my previous paper 112111 . 
I also introduced the notion of a regular violation of the fu- 
sion rules, if the fusion exponents ^„ p satisfy the inequality 
< ^n.p < which is necessary and sufficient to ensure that 
Ai and A2 are not governed by negative exponents, since, by 
definition, F„ should vanish in the limits Ai — > or A2 — >^ 0. 
For the cases p = I and p = n — I the fusion exponents are 
given by <^„,i = ^2 and <^„.„_i = For the case /:> = 1 it is 
assumed that one of the end points of the small velocity dif- 
ference is attached to one of the big velocity differences, that 
is Xi — X2. This assumption holds when working with par- 
tially fused correlations. It also holds when integrating F2 to 
compute the energy spectrum. 
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For the inverse energy cascade of two-dimensional turbu- 
lence, as was shown in my previous paper (iT], the fusion 
exponents are given instead by £,„,p = — C,n-p for 2 < p < 
n — 2. For the cases p = I and p = n — 1 the fusion expo- 
nents are given by ^„.„_i = ^„ — ^2 and ^„ 1 = 0. Again, the 
case p = n—1 requires the assumption that the large velocity 
difference is attached to one of the small velocity differences, 
that is x„ — x„_i. Note that, compared to the direct cascade, 
the roles of small and large velocity differences are switched. 

The fusion rules can be shown to be a consequence of the 
universality hypothesis, which we have introduced previously 
ll2Tll in a generalized form that is applicable to both the upscale 
and the downscale ranges of two-dimensional turbulence. Ac- 
cording to the universality hypothesis, the conditional gener- 
alized structure functions <!>„, given by 



<i>„({x}„,{Y}„„{w,r,Li,f) 



k=l 



w(xi,x'i,f) = W(i) 



(22) 



will honor the same symmetries, with respect to the point- 
pairs {X}„, as the unconditional correlations F„, in the asymp- 
totic limit where ||{Y}„,|| is situated between ||{X}„|| and 
the forcing scale £o- These symmetries are: incremental ho- 
mogeneity, incremental isotropy, and self-similarity with the 
same scaling exponents i^„. The ensemble here is constrained 
by a control condition imposed on the velocity differences 
w(xi,x'i.,f). 

It can be shown that the values of the fusion rules scaling 
exponents ^ can be derived as a consequence of this hypoth- 
esis of universality for both inertial ranges tf2Tl l26ll28ll . From 
the fusion rules we can then show that the integrals in the non- 
linear interactions term 0„f„+i are local, which implies that 
the scaling exponent of 0„F„+i is ^„+i — 1. If the fusion rules 
fail, then as long as the violation is regular, UV localiy is still 
maintained. Let us assume that p = ^p + A£,„,p in the down- 
scale range and ^„ p — ^„ — i^„_p + A^„ p in the upscale range, 
with A£,n p the corresponding perturbation. Then, IR locality 
locality is also maintained if the deviations in the fusion ex- 
ponents satisfy A(^„+i,2 + A^„+i_„_i > (downscale range) or 
A^n+1,2 + A^„+i_„_i < (upscale range), for n > 1. If both 
UV and IR locality are maintained, then the scaling exponent 
of 0„F„+i is still C„+i - 1- 

Knowing the scaling of the nonlinear term, it becomes pos- 
sible to compare it against the other terms of the balance equa- 
tions. As we have shown in our previous paper ll2Tll . if we as- 
sume that the forcing term is random-gaussian, then we may 
calculate the scaling exponent of Q,, and compare it against 
On^Vi+i- Thus, one can effectively circumscribe the boundary 
between the region 3n and the forcing range. In the present pa- 
per we will show how to compare the nonlinear term against 
the dissipation terms vJ„ and j3H„. In doing so, we deter- 
mine the position of the corresponding dissipation scales and 
may then show that they are positioned consistently with the 
anomalous sinks hypothesis. Thus, we will show that the 
anomalous sinks hypothesis can be derived as a consequence 
of the fusion rules hypothesis, which in turn follows from the 
proposed universality hypothesis. 



III. DISSIPATION SCALES FROM THE FUSION RULES 

The obvious problem with estimating the position of the 
dissipation scales by dimensional analysis is that by doing so 
we presuppose the validity of a correspoding anomalous sink 
hypothesis. The argument that we will use instead is similar 
to the argument of L'vov-Procaccia |27] for the case of the 
energy cascade of three-dimensional turbulence. We now re- 
view this method, show how it can be extended for the inverse 
energy cascade, and present an additional consideration that 
is very relevant in two-dimensional turbulence. 

Consider the case where in the correlation /^,'''(r,/?) be- 
tween n velocity differences, one of them is evaluated at sep- 
aration r and all others at R such that r <C <C ^q. Here, 
we assume that both r and R are still in the inertial range. 
For sufficiently small r, universal scaling fails as one leaves 
the inertial range and enters the dissipation range. The scale 

(n) 

£i,v where this cross-over occurs is dependent on R. Thus, 

one may define a dissipative length scale function r — il") (R) 
whose graph traces out the dissipative boundary of the en- 
strophy inertial range in the {r,R) plane. Similarly, in the 
case where £0 <^ r <^ R, one may define a large-scale dissipa- 
tion length scale function R = £^"^ (r) that marks the crossover 
point where the function fi" ^\r,R) enters the large-scale 

dissipation range when R exceeds (-["\r). The function ^|"''(r) 
sketches out the shape of the inverse energy cascade in the 
{r,R) plane. It should be remembered that a complete repre- 
sentation of the extent of one of the inertial ranges requires at 
least an 2n dimensional region 3„ C R^" for the correlation F„. 
Nevertheless, this two-dimensional representation is sufficient 
for our purposes, and certainly a step forward from the usual 
one-dimensional representation. 

In the enstrophy inertial range, when r is still in the inertial 
range, the function Fn^\r,R) evaluates according to the p=l 
fusion rule, and it is given by 

(23) 

Since we are ultimately interested in the dissipation length 
scales of structure functions, we may safely assume that one 
end point of the small velocity difference coincides with one 
of the end points of another velocity difference. The validity 
of the p ~ 1 fusion rule, used above, is contigent on mak- 
ing this assumption. When r enters the dissipation range, the 
dominant balance in the balance equations is 

vVl'^Fy'>ir,R)r 



R 



(24) 



and solving for F^^\r,R) yields 



F,i'Hr,R)^ 



vR 



(25) 



The dissipation length scale function is determined by match- 
ing the two asymptotic expressions for Fy\r,R): 



FliR) 



vR 



(26) 
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and solving for £i"){R). The next step is to determine the 

(n) 

length scale A„v where a cross-over to the dissipation range 
occurs when all velocity differences are shrinked simultane- 
ously. We estimate the location of this cross-over scale by 

solving the equation -^["'(Ai"') — Xu") with respect to A,!"'. If 
the enstrophy cascade forms successfully, then the dissipation 
length scale of the energy spectrum will be approximately lo- 

(2) (2) 

cated at Xin- , thus A„v is effectively an estimate of the Kol- 
mogorov microscale. Similarly, the dissipation scales of the 

standard structure functions S„ will be approximately at Xi") . 
In this sense, we say that Xjiy^ are the observable dissipation 
length scales. 

Although the observable dissipation scales are the ones 
that we are ultimately interested in, the shape of the function 

iuv{R) is also significant in the following sense: Suppose that 
all the velocity differences have been shrinked down to the 

scale Xu"), and we begin to stretch simultaneously n — 1 veloc- 
ity differences by a factor a > 1 to length aXi") while adjusting 

the remaining velocity difference to length £l^J{aXu"'^) such 
that we remain on the boundary between the inertial range 

(n) 

and the dissipation range. Then, the shape of the function £„y 
has to be such that we may remain on the boundary without 
being forced to increase the separation of the small velocity 
difference more than we have increased the other n — 1 ve- 
locity differences. In other words, we propose that in a valid 
inertial range it is necessary that 



(27) 



When this condition fails, we can have the inconsistent sit- 
uation of exiting the inertial range and entering the dissipa- 
tion range if we increase all velocity differences separations 
simultaneously by the same factor! This condition is obvi- 

(n) 

ously satisfied when ila- is a decreasing function. However, 

(n) 

it can also be satisfied when £uv is increasing, as long as the 

curve r = (R) remains underneath the line r = R. It will be 
shown later that this condition fails for the enstrophy range in 
certain cases. 

The same technique can be applied to the inverse energy 
inertial range. The evaluations of the function pj," ^\r,R) in 
the inertial range and the dissipation range are given by 

F2{R)F„{r) 



F!;"-'\r,R) 



Ft'\r,R)^ 



F2{r) ' 
R-^"'F„+i{r) 



(28) 



and the dissipation scale function is found by solving the 
equation 



(29) 



F2{4:\r))F„ir) (4"H0)-^"'^«+l(^) 
F2{r) ^ pr 

After solving the equation (["\xl"^) = Xl"\ we demand that 

(n) 

the shape of the function should satisfy the condition 

aA,W < 4")(«A,W) ,V« e 1). (30) 



We would Uke to stress again the significance of the form 

(n) (n) 

of the functions and i-^ . Consider for example the case 
of the enstrophy cascade. Suppose that we plot the lines r — R 

and r — £ul) (R) on the {r,R) plane. These lines intersect at the 

observable dissipation scale A^"' and they may or may not in- 

(n) 

tersect elsewhere. The region between the line r — £\J (R) and 
the line r — R can be thought of as a two-dimensional repre- 
sentation of the region where dissipation is negligible. Every 
point within that region represents a set of velocity differences 
configurations for which the correlations maintain self-similar 
universal scaling. To actually have an inertial range, the curve 
r = £^uv [R) has to be underneath the line r ^R from the point 
of intersection at the observable dissipation scale A]"' at least 
up until the forcing scale £q; the region between them needs to 
be reasonably inflated to allow local interactions to take place. 
A similar argument is applicable for the case of the inverse en- 
ergy cascade. 



IV. THE CASE OF THE ENSTROPHY RANGE 

We now employ the method described in the previous sec- 
tion to study the geometry of the enstrophy range in the (r,/?) 
plane. In this section, we let r\,„ represent the amount of the 
enstrophy that flows downscale. The value of Tj,,,, is the re- 
sponse of the system to the combined forcing and large-scale 
dissipation that drive the enstrophy cascade. A Reynolds num- 
ber can be defined for the enstrophy range as 



■J<-uv — . 



(31) 



We will consider two possibilities: an enstrophy inertial range 
with intermittency corrections and a logarithmic enstrophy 
range without intermittency corrections. Although the pos- 
sibility of intermittency corrections was recently ruled out by 
Eyink [43], it is still instructive to consider the possibility on a 
hypothetical basis. The reason for doing so is because we will 
show that in the case K" 1, corresponding to standard molec- 
ular diffusion, the geometry of the boundary between the in- 
ertial range and the dissipation range offers sufficient grounds 
to reject the intermittency scenario! Furthermore, we show 
that the logarithmic correction is needed to avoid contamina- 
tion of the inertial range from the dissipation range. In the 
case K" > 1 we show that any hypothetical intermittency cor- 
rections have to satisfy the constraint ^2 = 2 — ^2 < 2(k' — 1) 
to prevent dissipative contamination, although, as was shown 
by Eyink |43J, we have an additional reason for ruling out the 
intermittency corrections scenario altogether. 



A. Enstrophy range with intermittency corrections 



As a point of departure, we assume the evaluation given by 
the 1 /8 law that F-i{r) ^ {r]uvr^)- A derivation of the 1 /8 law 
is given by Lindborg [44] and Bernard 14511 . For the other 
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correlations we will allow the general form 



(32) 



where the scaling exponents ^„ satisfy ^2 —2 + 82 and ^„ = 
n — 5„ with 5,, > for n>2. It follows that the evaluation of 
the correlation F2^\r,R), when r is in the inertial range and 
the dissipation range correspondingly, reads 



F^'\r,R)^F2{r)^{r/£o)^^--\ri!J\)\ 



F^'\r,R) 



Fm 



(33) 



vR 



These evaluations coincide at the boundary between the iner- 
tial range and the dissipation range. For a fixed R, the cor- 
responding dissipation length scale £uv{R) is found from the 
following matching condition 



and solving for iuv{R) we obtain 



4 



(34) 



(35) 



We may transform this relation into an equation by introduc- 
ing a constant of proportionality via a critical Reynolds num- 



(77) 

ber „\, and define 



,{R) as 



iurjR) 



X 



in) 



(r?) 



l/(?2-2>c) 



(36) 



Recall that this length scale marks the point r = tuv{R) where 
one enters the dissipation range in a correlation of two velocity 
differences where one of them is held constant at R while the 
other is shrunk down to r. Implicit in our use of the fusion 
rules is the assumption that the two velocity differences have 
one point in common. It is remarkable that the dissipation 
scale li,v{R) is an anomalous function (i.e.. R dependent) even 
for F2. This should be contrasted with the energy range of 
three-dimensional turbulence where the dissipation scale is an 
anomalous function only for the higher order correlations F„ 
withn>2. iH. 

In the hyperdiffusion case K" > 1, we have 11,2— 2k < 
0, when the hypothetical intermittency corrections respect 
Eyink's constraint ^2 < 1 1/3 |46]. It follows that with increas- 
ing Reynolds number the separation of scales increases with 
the dissipation scale iuv{R) approaching 0. When fc < 1, with 
increasing Reynolds number, the dissipation scale iuv{R) di- 
verges to large scales instead of small scales, thereby demon- 
strating that such terms cannot provide a dissipation sink at 
small scales. The case K" = 1 with ^2 > 2 is interesting because 
it does not provide a well-behaved function either, and see be- 
low for more comments. Finally, when K = I and ^2 — 2, 
this evaluation is not valid, unless the logarithmic correction 
is taken into account. 



It should be remembered that the dissipation scale imiR) 
is not observed in the energy spectrum. The "standard" dissi- 
pation scale that we do observe can be found by solving the 
equation £uv{Xuv) = Kiv This leads to. 



in) 



in) 

2, Ml' 



l/(?2-2)c) 



where the scaling exponent A is given by 

2 C2-2()C+1) 82-2K 



1 



C2-2if C2-2if &2-2{K-lY 

and the solution is 

l/(&-2(/c+l)) 



in) 



(37) 



(38) 



(39) 



When ff > 1, the separation of scales between this dissipation 
scale A„v and the integral length scale £0 still increases with 
increasing Reynolds number, as we expect it to. For the case 
K" = 1 and ^2 — 2, this evaluation is identical to the dissipation 
scale evaluation that can be obtained by dimensional analy- 
sis. However, as we have noted in the previous paragraph, the 
function £uv{R), from which this evaluation has been obtained, 
is not valid in the absense of a logarithmic correction. 

The dissipation scale function £uv{R) is admissible if it sat- 
isfies the condition 



The condition holds if and only if 



aXuv 

^^-l+2/(f2-2;c) 



< 1, Vfl > 1, 



which is true if and only if 



(40) 

(41) 
(42) 

(43) 



Note that the same condition guarantees that goes to zero 
when the Reynolds number is taken to infinity. For fC > 1, 
any hypothetical intermittency corrections have to satisfy the 
inequahty ^2 < 2(k' — 1). 

The remarkable result is that in the physically relevant case 
K" = 1, the constraint imposed on any hypothetical intermit- 
tency corrections is ^2 > 2 (steeper than k^^). Because it 
contradicts the constraint ^2 < 11/3 of Eyink ll46ll . we con- 
clude that we may not have an enstrophy inertial range with 
intermittency corrections when fC = 1. In exact terms, we 
have shown that an enstrophy cascade with intermittency cor- 
rections will be destabilized if it is dissipated by standard 
molecular diffusion. It is an interesting coincidence, and one 
that may warrant some further reflection, that the slope l<r^ 
also occurs in the Tran-Bowman theory fl3l - [l5ll at the down- 
scale range, where the two-dimensional Navier-Stokes equa- 
tion does not have an infrared sink. 
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Although intermittency corrections do not seem to be for- 
bidden in the hyperdiffusion case, as far as stability with re- 
spect to dissipation is concerned, a more careful study [43] has 
akeady revealed additional constraints that exclude the inter- 
mittency scenario altogether, for the case of an enstrophy cas- 
cade with asymptotically constant enstrophy flux. Numerical 
simulations of the enstrophy cascade also corroborate the ab- 
sense of intermittency corrections (4^-^^. Nevertheless, this 
analysis sheds some light into the fundamental differences be- 
tween an enstrophy range dissipated by ordinary diffusion and 
one dissipated by hyperdiffusion. In the hyperdiffusion case, 
an enstrophy range that has already been somewhat destabi- 
lized into a steeper slope by the presense of a downscale en- 
ergy flux [21] will not be further disturbed by the dissipation 
range. In the ordinary diffusion case, we may find it more dif- 
ficult to obtain a robust enstrophy cascade with k^^ scaling, 
because any steepening deviation from that scaling caused by 
perturbation from the forcing term Qn or the sweeping term /„ 
would provoke further disturbance from the dissipation range. 
This is why it is very significant to investigate numerically 
whether an enstrophy cascade can exist under regular diffu- 
sion. 

To put these claims on solid ground, we must also consider 
the dissipation scale functions for the higher order correla- 
tions Fn. As usual, one velocity difference is shrinked down 
to scale r while the others are kept constant at scale R. It is 
further assumed that the small velocity difference is sharing 
an endpoint with at least one of the other velocity differences. 
Let us assume that the boundary between the enstrophy range 

and the dissipation range is located on the line r = i"uv {R)- To 

(n) 

find the function £„,, , we start from the matching condition 



F2iiln}mF,m {&{R)f''F„+,{R) 



(44) 



F2{R) vR 

in) 

and solve for {R). A simple calculation, given in detail in 
Appendix lAl gives: 

(n) 1 i/(?2-2<c) 

(45) 



t\R) 



4 



R 



31 



where the scaling exponents Xn are given by 

- C" + C2- 1 



— 



(46) 



Note that the use of matching conditions, corresponding to 
different values of n, may be introducing a different constant 
of proportionality, consequently the critical Reynolds num- 

bers 3?,,,t,v will depend on n. These higher order dissipation 
scale functions have the same dependence on the Reynolds 
number as £uv{R)- So, in that respect they behave similarly. 
The only difference is that the anomalous correction now 
scales with x„. 

The observable dissipation scales for F„ are found by solv- 
ing the equation £i"J{Xu")) = leading to 



(']) 



(47) 



Details are given in Appendix lAl Note that as the Reynolds 
number goes to infinity, the dissipation scale Xi") goes to zero, 
as it is supposed to, even for the case fc = 1 . Still, that does not 
necessarily mean that all these scenarios are admissible. The 
inertial range of the higher order correlations is admissible if 
it satisfies. 



(«) 



< 1, Va > 1. 



(48) 



This happens if and only if {x„ — l)lna<0, Vfl>l, and that 
requires 1 — x„ > 0. Note that ^„+i — ^„ > 0, by the Holder 
inequalities for a downscale cascade [21], which implies that 
(1 -x„)(C2-2ff) = 1 -2(f- (C„+i - C«) < 1 -2)f < 0,Vrc > 
1. It is therefore necessary that ^2 ^ 2k" < 0. For the case K > 
1, this condition leads to the same constraint ^2 < 2{k— 1 ) for 
the hypothetical intermittency corrections, which was derived 
from the admissibility condition on iuv{R). For the case, fc = 1 
the constraint 1^2 — 2k: < can never be satisfied when ^2 > 2, 
and that again excludes intermittency corrections. 



B. Enstrophy range without intermittency corrections 

We now repeat the previous analysis for the case C,„ — n and 
K" = 1 . Note that when K" > 1 , the evaluation of the dissipation 
scales in the intermittency case also applies to this case with 
subleading corrections introduced by the logarithmic correc- 
tion to the structure functions F„ which, to first approximation, 
we can safely ignore. When jc = 1, on the other hand, the role 
of these logarithmic corrections is crucial and has to be taken 
into account. 

We use the following evaluation for the correlations 



F„{R)-{rili'Rr[WoiR)r", 



(49) 



where a„ are the scaling exponents of the logarithmic correc- 
tion. The prediction of Falkovich and Lebedev ifsollsill is that 
a„ —2n/3>. To first approximation, we have disregarded the 
Bowman correction [52] to the logarithmic factor, which is 
negligible when R is not close to Iq, in terms of order of mag- 
nitude. The balance condition for the dissipation length scale 
of F2 leads to 



{r]lhuv{R)?[Hk,lluv{RW ^ 4v(^)'^[ln(4//?)]. 

VK 



(50) 



Note that without the logarithmic correction, the dissipation 
scale function £„v{R) cancels out completely from the match- 
ing condition. This means that the condition is always satis- 
fied, and therefore a hypothetical k^^ spectrum without log- 
arithmic correction would have to be part of the dissipation 
range for all wavenumbers k. The presence of the logarithmic 
correction however gives a solution for the dissipation scale 



,v{R)- Introducing a critical Reynolds number X 



2 ;;,„ we may 
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therefore define luv{R), as a leading order estimate, as: 



-- exp 



3/2 



3i! 



2,HV 



3/2 



(51) 



In the limit of infinite Reynolds number, the dissipation scale 
£uv{R) goes to zero, thus it is well-behaved. Solving the equa- 
tion £uv{^uv) = ^uv, the observable dissipation scale is found 
to satisfy 



3 r 



In 



1-1/2 



3/2 



{ri) 



(52) 



This equation cannot be solved in closed form, however when 
the separation of scales is large we get 



(53) 



consistent with the evaluation obtained from dimensional 
analysis. 

In general, the dissipation scale functions (ul) (R) of the cor- 
relations Fn can be derived directly from the matching condi- 
tion 




FiiR) 



vR 



(54) 



Solving for 4v'(^) we find, for the case of generalized, as 
opposed to Falkovich-Lebedev scaling. 



&iR) 



exp 



3^ 



l/«2 



51. 



2/fl2 



ln(^ 
R 



(55) 



where the logarithmic scaling exponents b„ are given by 

On+l — fl„+fl2 



fl2 



(56) 



The details are given in appendix|B] We stress again that this 
result is confined only to the special case fc = 1 in which the 
logarithmic factors are dominant, as discussed previously. For 
Falkovich-Lebedev scaling ifsol IsTll . the scaling exponent a„ 
reads a„=2n/3, which gives b„ = 3/2 for all n. For n — 2 and 

a„ = 2n/3, the equation for itv{R) reduces to our previous 
equation for luv{R)- 

The observable dissipation scale Xl,") is the solution of the 

equation l^^v (■^iv ' ) = X}") which simplifies to the following 
transcendental equation 





(57) 



and in the limit of large Reynolds numbers it vanishes as 

We can now confirm that the dissipation scale function 

£uv(R) satisfies the admissibility condition £l"}{aXi"}) < 

aX!,"} for all aX!'") £ {Xi") ,io). The argument is somewhat 
tedious, but given in detail in appendix |B] We find that the 

admissibility condition is satisfied if Xt') < ioexp{—b„a2/2). 
This condition is readily satisfied with sufficient separation 

between A„v and £(,. For Falkovich-Lebedev scaling, we 
have 02 — 2/3 and b„ — 3/2 and the condition reduces to 

x!,"^ < £o/ \f^- It should be noted that the proof in appendix 

|B]uses the claim that £iv^(fo) — £o which requires in turn that 
b„ > 0. In general, utilizing only mathematical, as opposed to 
physical considerations, it can be shown that the assumption 
that the logarithmic scaling exponents 02 and satisfy 02 >0 
and 03 > implies that b„ > for all n e N with « > 1 (see 
appendixicl). 

We conclude that in the physical case (f = 1, Kraichnan 
scaling is admissible but it is also necessary to have a loga- 
rithmic correction to allow the inertial range to form. In effect, 
the role of the logarithmic correction is to "inflate" the region 
between the curve r — (["J (R) and the line r = R. In the hy- 
pothetical case where the logarithmic correction is absent, the 
two curves will coincide and that won't leave any room for 
local interactions. 

It is interesting to note that when we have hyperdiffusion, 
the logarithmic correction is not needed to inflate the inertial 
range region. This points to an interesting difference between 
the case of physical diffusion and hyperdiffusion. In the hy- 
perdiffusion case, the inertial range will be reasonably inflated 
as soon as the energy spectrum has nearly converged to k^^ 
scaling, and we may expect that it will exhibit inertial range 
behavior before it converges completely. In the case of phys- 
ical diffusion, on the other hand, the enstrophy range will not 
begin to inflate until it converges to k^^ and begins converg- 
ing toward the logarithmic correction. In fact, it is possible 
that the absence of a region in the {r,R) plane in which local 
interactions can occur, when the spectrum is still steep, might 
make it impossible for the enstrophy range to converge at all. 



V. THE CASE OF THE INVERSE ENERGY CASCADE 

In this section, we let e„ represent the upscale energy flux. 
We begin with the evaluation given by the 3/2 law |45, sj, 
I54I1 . that Fj, (r) ~ e„ r, and for the other correlations we use 

Fn{r)^{e,rrrl\r/£^)'^"-"l\ (58) 
The dissipation scale is obtained from the matching condition 

F2{£^:^{r))F„{r) (r) 
F2{r) ^ J3r 

and that leads for n = 2 to 



(59) 



{eir£ir{r)f'\£ir{r)/l,)^^-^" ^ , (60) 
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It follows that £ir{r) is constant with respect to r with 
£ir, and the dissipation scale £{,- reads 



£ir 



where "Rj^. is the Reynolds number corresponding to the in- 
verse energy cascade given by 



l/(?2+2m) 



(61) 



l/3„-2/3+2m 



(62) 



This dissipation scale corresponds to the partially fused cor- 
relation F2 where one of the two velocity differences is being 
stretched towards larger scales while the other velocity differ- 
ence is held fixed at r. Note that, unlike the case of the enstro- 
phy cascade, here the dissipation scale £ir{r) is independent 
of r. Consequently, for all r in the inertial range, ^,v(r) coin- 
cides with the observable dissipation scale A,, , corresponding 
to the cross-over to the dissipation range when both veloc- 
ity differences are stretched simultaneously. We have a sim- 
ilar situation in the energy cascade of three-dimensional tur- 
bulence ll27ll . In the limit of the Reynolds number going to 
infinity, the dissipation scale also goes to infinity because 
^2 + 2m > 0. It is interesting to note that this also holds when 
m = 0. In other words, Ekman damping can provide a sink for 
the inverse energy cascade, whereas we've seen earlier that it 
cannot provide one for the enstrophy cascade. 

The dissipation scales that correspond to the correlations F„ 
with n > 2, on the other hand, are functions of r. Again, we 
consider the case where the velocity differences are partially 
fused, and one velocity difference is being stretched to large 
scales while all others remain constant at r, and we write the 
corresponding dissipation scale as £\"\r) = (-ir{r I ioY" . From 
this, we may rewrite ^2(4r'('")) 



2m 



-(r/io) 



(63) 



j3r(4"V))2'" 

Substituting that into the matching condition ( |59l ) reads 



Fn+l(r) 



F,{r) 



(64) 



After all the cancellations, we obtain the following constraint 
for the scaling exponents 

C„+i - (« + 1 )/3 = (C„ - «/3) - (C2 - 2/3) +x„(C2 + 2m). 

(65) 



Solving for x„ yields 

^ _ Cn+l ~ Cn + C2 ~ 1 

C2 + 2w 



(66) 



Note that ^,,+ 1 — Cw forms an increasing sequence, by the 
Holder inequalities 1I21I1 . and since X2 = 0, it follows that jc„ 
is an increasing sequence and therefore 



Xn > 0, Vn > 2. 



(67) 



Next, we solve the equation i["\xj'"'^) = Xj"\ and obtain the 
observable dissipation scales, which read 



Xir 



The dissipation scale function itself is admissible if it satisfies 
the condition 

^W^W) = a^"£^\xt^) = «^A(«' > aXl"\ V« < 1, (69) 
which holds if and only if 

(l-x„)lnfl<0, Vfl<l, (70) 



and this requires in turn that 

(2m+l)-(C„+i-C«) 



1 Xj-i 



C,2 + 2m 



> 0, Vn > 2. 



(71) 



This condition will be satisfied as long as the scaling expo- 
nents satisfy the following constraint 



C«+i ^ Cn < 2m + 1 , V« > 2. 



(72) 



From the Holder inequalities, we know that i^„+i — is 
either constant (for the case of no intermittency corrections) 
or increasing (otherwise) [21]. Therefore, in the presense 
of hypothetical intermittency corrections to the scaling expo- 
nents we should entertain the possibility that this inequal- 
ity could be violated for large enough n. If that were to hap- 
pen, then there would be no possibility of an inertial range for 
generalized structure functions F„ for large values of «, which 
would make for an interesting situation indeed. Here's what 
we can say with certainty: To violate this inequality for the 
least favorable case m = 0, we require multifractal contribu- 
tions to the inverse energy cascade with /z > 1 . In our previous 
paper ll2lll . we have shown that such contributions do not nec- 
essarily violate locality in the 0„F„+\ term. Furthermore, they 
do not violate the stability constraint h> 1/3 with respect to 
gaussian forcing either! If grows linearly as n — !> +00, then 
there exists an m such that the inequality can hold for all n. 
If, however, the growth rate of is faster than linear, then 
for any choice of m there is a no for which the inequality is 
violated for all n > «o- Consequently, as long as the question 
of whether the inverse energy cascade has intermittency cor- 
rections to the scaling exponents with superlinear growth 
rate remains open, the question of whether the inequality ( f72] i 
is broken at large n also remains open. 

Obviously, a possible violation of the inequality (f72] l. were 
it to occur, can be ameliorated by increasing m. On the other 
hand, it is reasonable to expect that a receding dissipation pro- 
file, that comes with increasing m, would eventually fail to 
dissipate the sweeping term effectively [16], leading to the 
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non-universal behaviour of Danilov-Gurarie |fTo| - [T2ll . Thus, 
the realization of a steady-state inverse energy cascade could 
require a balancing act in setting up the dissipation sink: a 
low-order m dissipates the sweeping term /„ more efficiently 
but makes it easier to break the inequality ( l72b . A higher-order 
m will delay the violation of ( |72] | considerably, but at the price 
of less efficient dissipation of the sweeping effect. 



VI. ON THE EXISTENCE OF ANOMALOUS SINKS 

The theoretical development presented in the last two sec- 
tions is based on two assumptions. First, on the fusion rules 
hypothesis, and second, on the evaluation that we have used 
for the generalized structure function F^. Naturally, both as- 
sumptions will eventually have to be justified. With respect 
to the first assumption, we will show in this section that if the 
fusion rules are not satisfied, then the dissipation scales will 
not be positioned in a manner that provides for anomalous 
sinks. Therefore, in cases where either the enstrophy cascade 
or the inverse energy cascade have been observed experimen- 
tally, the fusion rules hypothesis is at least as reasonable a hy- 
pothesis as the classical dimensional analysis arguments that 
have been put forth as a simpler explanation. 

Before we proceed with the main argument, we have to first 
take care of the following difficulty: In the argument of the 
preceding sections we have made the casual assumption that 
we may evaluate the generalized structure function ^3 using 
the 1 /8-law for the enstrophy cascade or the 3/2-law for the 
inverse energy cascade. The immediate concern is that these 
laws are usually established by assuming in advance that there 
is an enstrophy sink at small scales and an energy sink at large 
scales and that both sinks are anomalous. If we choose to 
make this assumption, then if we subsequently were to lever- 
age the argument of the preceding sections to prove that the 
sinks are anomalous, we would be engaging in a circular argu- 
ment. Fortunately, this does not affect an argument by contra- 
diction where we deliberately assume the existence of anoma- 
lous sinks and the violation of the fusion rules to deduce a 
contradiction. On the other hand, we would like to have a the- 
ory where it is possible to leverage the previous evaluation of 
the dissipation scales to investigate the existence of anoma- 
lous sinks. We will now explain how this can be done. 

Let us begin by considering, with no loss of generality, the 
case of the enstrophy cascade. The validity of the 1/8-law 
within an appropriate interval of length scales requires the fol- 
lowing three conditions: First, most of the enstrophy must 
flow downscale and most of the energy must flow upscale. 
Second, there must be a corresponding inertial range situated 
within a wide separation of scales between the forcing scale 
and the dissipation scale. Third, the positioning of the dis- 
sipation scale must be governed by a law such that for fixed 
viscosity it can dissipate any arbitrary amount of energy and 
enstrophy input. A similar set of conditions are needed for 
the validity of the 3/2 law, for the case of the inverse energy 
cascade. 

One possible interpretation of the argument of the previous 
two sections is that it has established that the first and sec- 



ond condition implies the third condition, under the assump- 
tion that inertial ranges, when they exist, satisfy the fusion 
rules. Let us consider this point carefully by summarizing 
the essense of that argument. First, we have assumed that the 
downscale range is an enstrophy cascade range, and the up- 
scale range is an inverse energy cascade range by presuming 
to know the scaUng of the generalized structure function Ft,. 
Thus, we have, so far assumed the first condition. We have 
also assumed the second condition when we asserted the exis- 
tence of an inertial range region governed by the fusion rules. 
However, contrary to a typical dimensional analysis argument 
which predetermines by hypothesis the location of the dissipa- 
tion scale, we did not also presume to know the extent of that 
inertial range. Instead, we have asked the theory to tell us if it 
is possible to have a region of scales in the (r,7?) plane where 
the dissipation terms of the balance equations are negligible 
relative to the local interactions term. 

We will now argue that the first condition can be established 
from first principles. As for the second condition, 1 will sug- 
gest that it is possible to deal with it provided that we adopt 
a more careful interpretation of the argument of the previous 
two sections. Then 1 will proceed with the anomalous sink 
argument itself. 

A. The logical structure of the argument 

The first condition, namely that most of the enstrophy must 
be constrained to flow downscale while most of the energy 
must be constrained to flow upscale, can be established from 
first principles, without assuming a priori the existence of 
anomalous sinks. The original argument, proffered as proof 
establishing this claim, was given a long time ago by Fj0rt0ft 
ifssll . It was later noticed |6, 56-58] that the claim underlying 
his argument, that the twin conservation of energy and enstro- 
phy is the sole determining reason that constrains the direction 
of fluxes in two-dimensional turbulence, is flawed. Instead, 
in a recent paper |59] co-authored with KK Tung, we have 
shown that it is the combined effect of the twin conservation 
laws and the mathematical structure of the dissipation opera- 
tor that constrains energy to go mostly upscale and enstrophy 
to go mostly downscale. 

Our key result, in that paper, is as follows: Let IlE{k) be 
the energy flux and let Ilg {k) be the enstrophy flux transfered 
from (0,^) to (^. +°o) per unit time by the nonlinear term in 
the Navier-Stokes equations. If the energy forcing spectrum 
FE{k) is confined to a narrow interval of wavenumbers ,A:2] 
such that FE{k) = when k ^ [^1,^2]^ then it can be shown 
for the forced-dissipative case, without making any ad hoc as- 
sumptions, that under stationarity, the fluxes II^ {k) and IIq [k) 
will satisfy the inequalities 

/ qYlE(q) dq <0,\/k> k2, (73) 
Jo 

q-^UGiq)>0,yk<ki. (74) 

The inequality ( |73] ) implies that the negative flux in the {Q,ki) 
interval is more intense than the positive flux in the {k2,+°°) 
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because the weighted average of TlE{k) gives more weight 
to the large wavenumbers. Thus, ( |73T l implies that energy 
fluxes upscale in the net. Similarly, (l74l ) implies that enstro- 
phy fluxes downscale in the net. It follows that we may justify 
to some degree the evaluations used for F% in the sense that 
they represent the only self- consistent scenario in which iner- 
tial ranges exist. 

Unfortunately, this still leaves open the question of the suc- 
cessful existence of inertial ranges where universal cascades 
are to be situated. In the argument of the previous two sec- 
tions, taken according to the more obvious interpretation, we 
seem to assume without proof the existence of a forcing- 
dissipation configuration that allows cascades to exist. How- 
ever, we can view the same argument in an entirely differ- 
ent context if we situate it within the broader argument, out- 
lined in greater detail in the conclusion of the present paper, 
of which it is a part. 

The structure of the broade r arg ument, in summary form, 
consists of 2 basic steps ifTsi Il9l I21II : First, we show that 
the balance equations have two homogeneous solutions, cor- 
responding to an energy cascade and an enstrophy cascade 
for both directions, and a particular solution driven by forcing 
and sweeping. The flux inequalities ( TTSl l and ( |74] l establish the 
downscale enstrophy cascade and the upscale energy cascade 
as the leading homogeneous solutions. Second, to show the 
conditions necessary for the existence of cascades with uni- 
versal scaling, we require the existence of a region C M^" 
where the homogeneous solutions dominate the particular so- 
lution in the generalized structure function F„{{%]„,t) for all 
{X}„ e An and for each choice of n e N with « > 1 . A sep- 
arate effect is the distortion of the homogeneous solutions by 
the dissipation term. This defines a separate region 23„ C M^" 
where this distortion effect is negligible. For the successful 
formation of a cascade we need to have a measurable region 
3„ = A„ n 23„ with non-zero measure wherein the regions A„ 
and 23„ overlap. This overlap region 3„ is a multidimensional 
representation of the extent of the inertial range for the corre- 
sponding cascade, whereby the leading homogeneous solution 
dominates the particular solution. 

The point is that both regions A„ and S„ are to be calcu- 
lated separately first, and the existence of the cascade itself 
is then decided by the existence of an overlap 0„ between the 
two regions A„ and 2?„. So, we can take the standpoint that 
the argument of the previous two sections is the determina- 
tion of the region A„, whereas the argument of my previous 
papers lfT6[ I21I1 concerns the determination of the region 23,,. 
From this standpoint, we can argue that in the previous two 
sections we were investigating the dissipative corrections to 
the homogeneous solutions independently of whether or not 
the homogeneous solutions actually dominate over the par- 
ticular solution. As we have shown previously |18, 19, 21], 
the scaling exponent ^3 evaluations as ^3 = 1 for an energy 
cascade and Cs — 3 for an enstrophy cascade are an inherent 
characteristic of the homogeneous solutions themselves that 
is unrelated entirely from the anomalous sink question. 



B. Anomalous sinks from the fusion rules 

With the above prefatory remarks under consideration, let 
us now proceed with the anomalous sink argument. Let E{k) 
be the energy spectrum, and suppose that a certain amount of 
energy and enstrophy is dissipated at large scales and another 
amount at small scales. The energy dissipation rates are given 
by the following integrals. 



e,v==2)3 / k-^"'E{k)dk, 
Jo 

£„,. = 2v/ k+^''E{k)dk, 
Jo 



(75) 



and under steady state they satisfy e = e„ + £,„,, where £ is the 
rate of energy input. Here, e„ is the energy dissipation rate 
at large scales and £„,, is the energy dissipation rate at small 
scales. Similarly, the enstrophy dissipation rates are given by. 



T7,> = 2J3 / k-^'"+^E(k)dk, 
Jo 

T],,,, = 2v/ k+^''+^E{k)dk, 
Jo 



(76) 



and, again, under steady state they satisfy Tj = 77,, + TJ,„., where 
Tj is the rate of enstrophy input. Here, Tj,> is the enstrophy 
dissipation rate at large scales and tj,„, is the enstrophy dissi- 
pation rate at small scales. 

We begin with the case of the inverse energy cascade. Since 
the effect of the sink at small scales can be safely ignored at 
large scales, the dominant contribution to the energy dissipa- 
tion rate is given by the integrals 



Sir = 2/3 



1/A„. 



k-^'"E{k)dk + 2P 



-2m 



E{k)dk 



2B / k-'-"'Eik)dk. 
1/4 



(77) 
(78) 



Here A,> is the observable dissipation scale at large scales. For 
the case of the inverse energy cascade we have found that 
C.ir{r) is independent of r with £,v(r) = A,>. We have also 
shown that 



3^, 



l/(?2+2m) 



-l/(?2+2m) 



(79) 



The key point is that the scaling exponent ^2 in the above 
equation originates from the order R velocity difference sepa- 
ration in the {n,p) = (2,1) fusion rule, which generally scales 
as /;?2-&,i xhe fusion rule hypothesis gives ^2.1 —0 leading 
to the above estimate of A,>. If this fusion rule were to be vio- 
lated, the scaling exponent ^2 — ^2.1 would replace ^2^ and the 
dissipation scale A,v would then follow 



-l/(?2-&.l+2m) 



(80) 



In the limit of extending the separation of scales in the inverse 
energy cascade, we have £0 0+ which kills the third inte- 
gral and j3 ^ which gives A,> +00 which kills the first 
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integral. It follows that the dominant contribution comes from 
the second integral. When we substitute E{k) ^ A:^''^?^) [q 
the second integral we get: 



/t-2'«;t-(l + ^2't//t- j3(l/A,v)-(?2+2'") (81) 



'l/l,> 

PX^2+2m ^ ^l-(?2+2m)/(?2-fcj+2m) 



(82) 



It follows that the anomalous sink hypothesis for the inverse 
energy cascade which requires that Sir be independent of /3 
is satisfied if and only if the fusion rules hypothesis ^2,1 =0 
holds. In other words, we cannot have an anomalous sink if 
the fusion rules hypothesis is violated. It is interesting to note 
that large m has a tendency to mitigate the impact on the sink 
anomaly of a small discrepency between ^2 — ^2,1 and 1^2 ■ That 
is, when m — > +0°, £,-,. becomes independent of j3 regardless 
of the value of the scaling exponent ^2.1- 

A similar argument can be employed to establish the exis- 
tence of anomalous enstrophy sink at small scales. The domi- 
nant contribution to the enstrophy flux is given by the integrals 

Tj„v = 2v/ k^''+^E{k)dk + 2v k^''+^E{k)dk (83) 
Jo J 1/(0 



-2v I k^''+^E(k)dk. 



(84) 



Here A„v is the observable dissipation scale at small scales. 
For the case of the enstrophy cascade with hyperdiffusion (i.e. 
K" > 1) we have found that £uv{R) is anomalous, in the sense 
that it is R dependent and reads 



3^: 



(ri) 
'2,iiv 



l/(S2-2lc) 



(85) 



As long as 1^2 — 2 k: < 0, the observable dissipation scale X^y is 
the solution of the equation £uv{Xuv) = Aj,,, and it reads: 



Xuv 



in) 



5i: 



2,HV 



l/(f2-2(ic+l)) 



,-l/(&-2(ic+l)) 



(86) 



Again, the scaling exponent ^2 now originates from the order 
r velocity difference separation in the {n,p) — (2,1) fusion 
rule, which scales as r''^ 1 , consequently, under violation of 
the fusion rules hypothesis, the scaling exponent ^2,1 would 
replace ^2, and the dependence of A„v on v would be given by 



,-l/(fe.l-2(^+l)) 



(87) 



Now, by repeating a similar line of argument as in the previous 
case of the inverse energy cascade, we argue that in the limit of 
extending the separation of scales in the enstrophy cascade the 
dominant contribution to rj^y comes from the second integral. 
Substituting the enstrophy spectrum E{k) ~ A;^'^^^^) to the 
second integral gives: 

TJ„v - V ^ v(l/A„,)2'^+2-& (88) 

- vX^r^'-"^^^ ~ ^,l-(?2-2(rc+l))/(^2.l-2(/c+l))^ (gg-) 



Consequently, we find again that the anomalous sink hypothe- 
sis for the enstrophy cascade which requires that rjuv be inde- 
pendent of V is satisfied if and only if the fusion rules hypoth- 
esis ^2,1 — C2 holds. Again we see that large K hyperdiffusion 
helps the enstrophy sink to remain approximately anomalous 
even under a small deviation of (^2,1 from ^2- 

The careful reader will note that the evaluation of A„v used 
in the argument above is not valid when K = 1 and ^2 = 2. 
We have already shown that in this case, Xuv is given by the 
transcendental equation 



In 



X„v 



1-1/2 



3/2 



2,HV 



for the case of Falkovich-Lebedev scahng 
that 



(90) 



IT. It follows 



(iA„ 



/-3/2r 



n(4/A„„)]'/^ 



(91) 



and thus, repeating the calculation of Tjj,,,, with the inclusion 
of the logarithmic correction to the energy spectrum E{k) ~ 
/fc-3[ln(Mo)]-i/^ 



now gives: 



k'^k-^[ln{k£o)r^^^ dk 

v(l/A„)2[ln(4/A„„)]-'/^ 
v[(l/A„,)3]2/3[ln(£o/A„v)]-'/' 



const . 



(92) 

(93) 

(94) 

(95) 
(96) 



which shows again that we have an anomalous sink. 

Incidentally, it is an interesting exercise to consider the case 
of arbitrary logarithmic scaling and derive the necessary and 
sufficient condition on the exponents fl„ for the existence of 
an anomalous energy sink. In general, the energy spectrum 
scales as E{k) k'^[ln{k£o)]"2'\ and it follows that 



1 / A-Hv 

nuv-^v I k'^k-^[\n{kloyf^-Uk 

^v(l/A„)2[ln(4/A„)]"^"'- 
The transcendental equation for Aj„. now reads 

1/^2 



(-0 



2/a2 



In ^ 



4 
A,, 



I -In 




(97) 
(98) 



(99) 



which gives 

A) ^V-1/2 
A,/v / 

It follows that 



In 



A,/v 



b2- 1)02/2 



r?„v - v{v-'/2[ln(4/A„,)](''2-')"2/2}2[in(^^/^^^,)]«2-i 
~ [\n{£^,/Kv)f^-'^"^+^"^-'^ ^ [ln(^o/A,„)]"2'2"'. 



(100) 



(101) 
(102) 
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Since 02^2 ^1 = 03 — 1, we see that an anomalous enstrophy 
sink exists if and only if 03 = 1 . This is obviously consistent 
with the Falkovich-Lebedev claim ifsollsil] that a„ = 2n/3. 

It is possible to combine this argument and the previous ar- 
gument to generalize this result for the case K" > 1, provided 
that the logarithmic corrections are taken into account in cal- 
culating the location of the dissipation scale A„v. This gives 
the enstrophy dissipation rate rjuv as: 

r?„v - v^-(?2-2('^+i))/(«2^i-2(^+i))[ln(£o/A„v)]"'"'. (103) 

Again, the fusion rules hypothesis rules out the power-law de- 
pendence of rjiiy on V. However, the condition 03 = 1 is still 
needed to rule out a logarithmic dependence of tj„v on v. Fu- 
ture work should thus take another good look at these loga- 
rithmic scaling exponents. In connection with this remark, it 
is worth noting that, in recent papers |60, 61], Tran claims that 
he proved that it is not possible for the enstrophy cascade to 
have an anomalous enstrophy sink. However, the range of ap- 
plicability of his result is confined only to the case of freely 
decaying two-dimensional turbulence, with K" = 1 and no sink 
at large scales (i.e. j5 = 0). As such, it is not applicable under 
two sinks or for the forced-dissipative case. 

Although the case of the energy cascade of three- 
dimensional turbulence is not the primary focus of this paper, 
it is easy to see that our argument can be repeated to derive 
the dependence 

e _ ^,l-(f2-2(ic+l))/(fei-2(;c+l))^ (jQ4^ 

between the energy dissipation rate e and the viscosity V. 
Under the regular violation of the {n,p) = (2, 1) fusion rule, 
where ^2,1 satisfies < ^2.1 < we see that the energy dis- 
sipation rate e will vanish in the limit v — > 0+, suggesting an 
insufficiently powerful energy sink. If and when such a viola- 
tion does occur, the energy spectrum E (k) might respond by 
reducing the value of the scaling exponent ^2 until it converges 
to the value of £,2,1- That would lead to a shallower slope for 
the energy spectrum ^(fc). Such a response would not violate 
locality |21] and it is therefore reasonable to expect that the 
energy would still trickle down via local interactions and dissi- 
pate when dissipation is strong enough, thus maintaining sink 
anomaly in the intermediate asymptotic sense. Our argument 
shows that a regular fusion rules violation combined with in- 
termediate asymptotic sink anomally must lead to a shallower 
energy spectrum. This is consistent with a bottleneck-type 
pile up of energy near the interface between the inertial range 
and the dissipation range. Thus, a question that should be 
considered by future work is whether the well-known bottle- 
neck problem ll62il63ll of three-dimensional turbulence can be 
understood in terms of a regular violation of the fusion rules. 

VII. CONCLUSION 

Let us now reconsider, in broad strokes, the logical struc- 
ture of the overall argument, as it pertains to the cascades of 
two-dimensional turbulence. The point of departure are the 
generahzed balance equations. First, we observe that these 



equations have a general solution that consists of two homo- 
geneous solutions that represent correspondingly an enstrophy 
cascade with 1^3 = 3 and an energy cascade with Cs = 1 uM, 
and a particular solution driven by forcing and sweeping. The 
determination of the ^3 exponents for the two homogeneous 
solutions does not use the anomalous sink hypothesis. It fol- 
lows instead from a self-consistency argument. Second, we 
derive the constraints ( |73] ) and ( l74l i on the fluxes TlE{k) and 
11(3 (A;), which imply that the only self-consistent manner in 
which these cascades can manifest themselves is as a down- 
scale enstrophy cascade and an upscale inverse energy cas- 
cade 1I59I1 . It is possible, for the case of finite viscosities that 
side by side with the dominant cascades there may be a sub- 
leading downscale energy cascade and a subleading upscale 
enstrophy cascade ifisll . However, a constraint on the permit- 
ted magnitude of the energy and enstrophy fluxes associated 
with these subleading cascades, prevents them from having an 
observable effect on the energy spectrum lfl9ll64ll . 

As was discussed in my previous papers i[il[il[ill2il], 

in order for either the inverse energy cascade or the enstro- 
phy cascade to exist, there must be a region A„ C M?" where 
the corresponding leading homogeneous solution dominates 
the particular solution driven by sweeping and forcing. This 
is the so-called stability condition. Furthermore, the dissipa- 
tion terms effectively act to distort the homogeneous solutions 
within a certain dissipative region. If !B„ C M?" is the region 
where such dissipative effects on the leading homogeneous 
solution are negligible, then our main requirement for the ex- 
istence of a cascade is that there should be a measurable over- 
lap dn = -AnD Tin between the two regions A„ and 23,, with 
non-zero measure. The region 3„ is thus a multidimensional 
representation of the extent of the inertial range associated 
with the generalized structure function ¥„ . 

We may now pose the following question: Does the struc- 
ture of the mathematical theory of turbulence allow the exis- 
tence of cascades with universal scaling? In other words, if we 
require the cascades to have universal scaling exponents, can 
the region 2in have a non-zero measure? The argument in re- 
sponse to the question now runs as follows: First, the demand 
that the cascades must have universal scaling implies that the 
cascades, if they exist, must satisfy the fusion rules hypothe- 
sis. As things stand now, this first step is the incomplete link 
in the overall chain, as far as two-dimensional turbulence is 
concerned. This is discussed further below. 

Second, from the fusion rules hypothesis, we establish the 
locality of the local interactions term 0„F„+i of the balance 
equations and combined with the assumption of random gaus- 
sian forcing we then derive the stability conditions of either 
cascade with respect to forcing perturbations |21]. The ques- 
tion of stability with respect to sweeping is still an open prob- 
lem [16]. Both kinds of stability are necessary conditions for 
the existence of cascades with universal scaling, because if 
either type of stability were to fail, then we would reach a 
contradiction with the assumption, in the beginning of the ar- 
gument, that the cascades have universal scaling. Such contra- 
diction would imply the absence of universal scaling. We have 
shown |21] that the inverse energy cascade is stable with re- 
spect to perturbations of the forcing statistics. The enstrophy 
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cascade is only marginally stable, with the stability contingent 
on the smallness of the downscale energy flux. 

Third, from the fusion rules hypothesis and locality (a con- 
sequence of the fusion rule hypothesis) we have the argument 
of the present paper which leads to establishing the location 
of the dissipation scales, and, consequently, the existence of 
anomalous sinks of enstrophy and energy at small and large 
scales correspondingly. Our argument, more generally, estab- 
lishes the relation between the energy dissipation rate e,,- and 
the enstrophy dissipation rate and the corresponding fric- 
tion parameters j3 and v in terms of the scaling exponents Qi 
and ^2,1 • As a matter of fact, we have seen that the scope of the 
argument is broader: we have derived the shape of the bound- 
ary separating the inertial range from the dissipation range in 
the {r,R) plane representation for every generalized structure 
function ¥„ for any n > 2. 

As it stands now, the above argument has two loose ends. 
First, we have not worked out the problem of stability with 
respect to the sweeping interactions. As we have pointed out 
previously LI61I . this problem, which was first recognized long 
ago by Kraichnan ll65ll . remains an unresolved issue that af- 
fects every theoretical approach to the problem of turbulence 
that has been proposed to date. Second, the first step of the 
argument, establishing the fusion rules hypothesis from the 
assumption of universality, has not been done in a completely 
satisfactory manner. 

This last point requires further clarification. The currently 
used formulation of the hypothesis of universality is that the 
fundamental symmetries of the inertial range, namely incre- 
mental homogeneity, incremental isotropy, and self-similarity, 
defined in terms of the ensemble-averaged generalized struc- 
ture functions, will survive imposing a symmetry-breaking re- 
striction on the ensemble averaging. Thus, the conditional 
generalized structure functions <!>„, defined via a conditional 
ensemble average of products of velocity differences incorpo- 
rating the restriction, will still have an inertial range where 
the same symmetries are satisfied, as long as the symmetry 
breaking restriction is introduced at a scale that lies between 
the forcing scale and the scale defined by the velocity differ- 
ences that are being averaged. Now, there is a choice: The re- 
striction can be placed on the forcing field or on the velocity- 
difference field. If the restriction is placed on the velocity- 
difference field, then the fusion rules hypothesis is easily de- 
rived from the universality hypothesis. This was noted previ- 
ously by L'vov and Procaccia ll26ll28ll in the context of three- 
dimensional turbulence, and I used the same argument to pre- 
dict the fusion rules for both cascades of two-dimensional tur- 
bulence ll2ll] . 

Now, this is hardly a surprise. From a physical stand- 
point, the hypothesis of universality, formulated in terms of 
a velocity-difference constraint on the ensemble, means that it 
makes no difference to the symmetries at a certain scale, if that 
scale is forced directly by gaussian forcing or indirectly by 
the velocity field at a larger (or smaller, for the inverse energy 
cascade) scale. In other words, we are already, in a certain 
vague sense, hypothesizing the existence of a local cascade. 
Thus, the locality and stability arguments, given previously 
ll2lll . and the arguments of the present paper, essentially es- 



tablish a required consistency between the existence of local 
cascades and the governing Navier-Stokes equations. This is 
useful, but not good enough. 

To establish the fusion rules hypothesis in a more definitive 
manner, we should start from a weaker universality hypoth- 
esis in which we merely require that the scaling laws should 
not be sensitive to forcing perturbations, and nothing more. 
This is, after all, what we really intend to say when formu- 
lating the universality hypothesis. Arguments of this nature 
have been used to establish the fusion rules hypothesis for the 
passive scalar problem |34] of Kraichnan |66] and the case 
p = 2 has been previously proved for the energy cascade of 
three-dimensional turbulence I23l - l25ll . This type of stronger 
argument is currently missing for the case of the cascades in 
two-dimensions. 

This is why the mathematical argument, both in this paper 
and in the previous paper fl)^, was carefully formulated to 
account for the case where the fusion scaling exponents ^ 
get to deviate from the values predicted by the fusion rules 
hypothesis. For example, in the previous paper ||2lll I inves- 
tigated the effect a violation of the fusion rules will have on 
the locality of nonlinear interactions in the inertial range, and 
derived the required inequalities on the fusion scaling expo- 
nents ^„ p for preserving locality. In the present paper, I de- 
rived the general dependence of the downscale enstrophy dis- 
sipation rate rjuv on the viscosity v and the upscale energy 
dissipation rate £,> on the viscosity j3 in terms of the scaling 
exponent ^2 and the fusion scaling exponent ^2.1- For the case 
where ^2,1 satisfies the fusion rules hypothesis and 03 = 1, rj^v 
and £,, become independent of the viscosities v and j3, which 
implies the existence of anomalous sinks. However, it should 
be stressed that our analysis was deliberately done in a general 
way to cover the case where the fusion rules hypothesis fails. 

With all the aforementioned concerns notwithstanding, I 
believe that a big picture is beginning to emerge, which can 
sharpen our physical intuition towards a deeper understand- 
ing of the cascades of two-dimensional turbulence. Already, 
from the present argument, we have gained some interesting 
insights about two-dimensional turbulence, which shall now 
be summarized. First, we have shown that it is impossible 
to have an inertial range in an enstrophy cascade dissipated 
by molecular dissipation (k ~ 1), if the enstrophy cascade 
has intermittency corrections, or if, for some other reason, it 
deviates from Kraichnan scaling into steeper slopes. Conse- 
quently, it may be easier for the enstrophy cascade to con- 
verge under hyperdiffusion as opposed to molecular dissipa- 
tion. Second, we have shown that the logarithmic correction 
plays an essential role in ensuring that the inertial range of 
the enstrophy cascade is not entirely contaminated by dissi- 
pation, when fc = 1 . Third, for the case of the inverse energy 
cascade, we have shown that if there are intermittency correc- 
tions to the scaling exponents , then the scaling exponents 
must satisfy the inequality ^„+i — ^„ < 2m + 1 , Vn > 2 , with 
m being the order of the hypodiffusion, in order for all gen- 
eralized structure functions F„ to have an inertial range. Last 
but not least, we have not only shown that the fusion rules 
imply the existence of anomalous sinks; we have also estab- 
lished that a possible small violation of the fusion rules can be 
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compensated for by increasing the orders K and m of hyper- 
diffusion and hypodiffusion correspondingly. For the case of 
three-dimensional turbulence we have also pointed to a possi- 
ble connection between the bottleneck problem ll62i l63ll and a 
possible temporary regular fusion rules violation. 

To conclude, it should be noted that our work has focused 
strictly on the mathematical problem of two-dimensional tur- 
bulence that is precisely two-dimensional. More realism 
would require the introduction of three-dimensional effects 
at small scales, and still unaddressed remains the problem of 
how these effects impact on the dynamics of two-dimensional 
turbulence. In connection with this very broad question, we 
can make the following remarks. First, the reason why we 
have a sink at large scales is precisely to account for Ekman 
damping, which is a three-dimensional large-scale frictional 
effect, known to exist for the case of atmospheric turbulence 
iIHtIi . a similar frictional effect exists in soap-film experi- 
ments of two-dimensional turbulence ll68ll . Second, at small 
enough scales, one expects a change in the dynamic of the 
nonlinear term. Based on what we have learned about local- 
ity from my previous paper ll2lll . I expect that the dynamic of 
the two-dimensional inertial ranges will be shielded from the 
nonlinear small-scale three-dimensional effect. However, the 
three-dimensional effect could change the nature of the small- 
scale dissipation. Further speculating on this issue is beyond 
the scope of this paper Let us just emphasize again that, based 
on our results here, we have learned that it makes a signifi- 
cant difference whether the enstrophy cascade is dissipated by 
diffusion or hyperdiffusion. As an anonymous referee noted, 
"the short-wave dissipation is commonly considered as a pa- 
rameterization which has to be designed so as to dissipate the 
enstrophy flux. In contrast, in the 3D case we believe that 
the true dissipation operator is known". The fact that we do 
not know the true dissipation operator in two-dimensional tur- 
bulence, in light of our results, is thus revealed to be a more 
urgent problem than is commonly realized. 



Acknowledgments 

The author wishes to thank an anonymous referee for his 
constructive remarks. 



Appendix A: Dissipation scales for an enstrophy cascade with 
intermittency corrections 

In this appendix we give the detailed derivation of the dis- 



sipation scales lul) (R) and Xi,"^ for the case of the enstrophy f 
cascade with (mostly hypothetical) intermittency corrections. 

To find flie function il"), we write el"){R) = 4v(^)(^/4)"", 
and proceed to evaluate a„ using the matching condition 
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The factor F2(C^(/?)) can be evaluated with judicious use of 
the properties of the dissipation scale functions as follows 



vR ^ ' ' 
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We substitute this into the matching condition and obtain 
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If we employ the general evaluation of F3 and F„ we note that 
the factors {£"^^,{R))'^'', vR, and rjul^ r cancel completely. We 
are left with the balance 



(A6) 



that implies the scaling exponent equations 

C„+i - (n + 1) = - n - (C2 - 2) + fl„(C2 - 2ff). (A7) 
Solving for a„ we get 

C«+i ^ C« + C2 ^ 3 



(A8) 



We may now go back and write the dissipation scale function 
eil){R) as follows 
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where the scaling exponents x,, are given by 

_ 2 Cn+l ~ C» + C2 ~ 3 

_ Ci+l ~ Cn + C2 ~ 1 



(AlO) 



The observable dissipation scales for F„ are found by solv- 
ing the equation — X^"}, leading to 
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and thus. 
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A simple calculation gives: 

(l-x„)(C2-2jf) 



(A13) 



= (C2-2-c) 



(C2-2>c)-(C„+i-C« + C2-i) 



1-2k:-(C„+i-C«), 



(A14) 
(A15) 



and it follows that the dissipation scale Ajl"^ is given by 

(r,) 1 l/[l-2)f-(?„+i-«] 
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Appendix B: Dissipation scales for an enstrophy cascade without intermittency corrections 

In this appendix we give the detailed derivation of the dissipation scales il^J (R) and x!,"^ for the case of the enstrophy cascade 
without intermittency corrections for the special case where K =1. Then, we confirm the admissibility condition. 

To find the function we start by substituting F„(R) ~ {rill^R)"[ln{£o/R)]"'' to the matching condition 



F2{iuXR))Fn{R) _ {C'iR)fF„+r{R) 
which gives 

(?7,l/'/?)2[ln(4//?)]«2 
After a few immediate cancellations we get the condition 
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from which we find that 
In 
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Solving for the dissipation scale iuv{R), and introducing a critical Reynolds number 3?2 «v' ^® proportionahty relation 

into an equation, gives the main result that reads: 
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with bn given by 

On — • 
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The observable dissipation scale Ai"^ is the solution of the equation tuv{^uv') = ■^mv\ which reads 
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and then simplifies to the following transcendental equation 
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We may now confirm that the dissipation scale function fi"' (R) satisfies the admissibility condition fi"' {aXu") ) < aXu") for 
all aXi,'^ G {X!i"} ,£q). First, we note that il!l){£Q) = £q (requires that b„ > 0) and £m}{X,n^) = A^"'. This means that the curve 
r = (-'uv{R) intersects the line r — R al the endpoints of the inertial range R ~ xl") and R = Iq. The admissibility condition 
requires that when R moves from Ai"' to Iq the curve r = il"^{R) must remain under the line r = R throughout the interval 
R e (Ativ^^o)- To confirm this behaviour, consider the partial derivative of £m){R) with respect to R: 
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All factors above are positive, including the logarithm 
ln{£()/R) (since R < £q), except for the last factor that can be 
either positive or negative. To satisfy the admissibility con- 
dition, it is sufficient to require that £^uv{R) be decreasing at 
R = Xi"^ . This occurs if and only if 




-7^ -^">0, 



(B13) 



which is satisfied if and only if X^^ < £oexp(— Z7„a2/2). For 



R > Xuv', the curve of £["J{R) will continue to decrease for a 
while, and will then turn around to increase and catch up with 

the line r = R, since £l^) {£o) = £q. 



Appendix C: Inequalities on logarithmic scaling exponents 

In this appendix we show that the assumption that the loga- 
rithmic scaling exponents 02 and 03 satisfy 02 >Q and aj, >Q 
implies that b„ > for all n G N with « > 1 . This is a conse- 
quence of the inequality fl„+i — a„ >a„— a„_i which can be 
established by an argument that is analogous to the one used 
to establish similar inequalities between the regular scaling 
exponents C« iSlliliilTQll. 

We begin by defining 5^{R) as the absolute value of the 
scalar vorticity difference: 



5CW = |C(x + /?e,r)-C(x,f)|, 



(CI) 



where x £ is given and e is a unit vector. The proof is 
based on the following two assumptions: (a) For a downscale 
cascade, in the limit £0 °°, ^C{R) scales as {[8^{R)]") ^ 



[ln{£Q/R)]"". (b) For finite £0 there is a range of scales where 
the above scaling law continues to hold as an intermediate 
asymptotic. 

Let p,q e (1,+°°) with l/p+ l/q = 1, and let (j),\jf 
be two random variables with (j) > and xj/ > 0. The 
Holder inequality for ensemble averages states that (0)/) < 
{0P)^^'' {-^1) For 79 = ^ = 1 /2 it reduces to the Schwarz 
inequality: {(j)yf)^ < We begin by choosing (j) = 

[SCl^)]'""''/^ and y/ = and employing the 

Schwarz inequality. It follows that 



= ([5cwr')([5cwr'). 
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and therefore 
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To satisfy this inequality under the limit £0 — °° we require 
2a„ — fl„_i — fl„^i < which is equivalent to a„^i — a„ > a„ — 

Let us assume now that 02 > and 03 > 0. It follows that 



b„ = 

^ a3- 02 + 02 _ ^ ^ Q 
~ fl2 ^2 ' 

which concludes the proof. 
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